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Abstract 

We compute numerically the full retarded Green's functions for conserved currents 
in the shear channel of a (2+l)-dimensional field theory at non-zero temperature and 
density. This theory is assumed to be holographically dual to a non-extremal, elec- 
tric Reissner-Nordstr0m AdS4 black hole with planar horizon. Using the holographic 
description we obtain results for arbitrary frequencies and momenta and survey the de- 
tailed structure of these correlators. In particular, we demonstrate the 'repulsion' and 
'clover-leaf crossing' of their poles and stress the importance of the residues at the poles 
beyond the hydrodynamic regime. As a consistency check, we show that our results 
agree precisely with existing literature for the appropriate quasinormal frequencies of 
the bulk theory. 
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1 Introduction 

The AdS/CFT correspondence is a valuable tool in the study of strongly-coupled systems. 
For example, the real-time correlators of such systems can be extracted from classical gravity 
PP, E] . By tuning the physics of the bulk theory, the boundary theory can be made to mimic 
the qualitative features of many condensed matter theories. Key examples are the superfluid 
phase transition [3J SI [5] and a Fermi surface [SI El El- See [9], [10] and [H] for reviews. 

Many previous studies have restricted to calculating just the locations of poles in re- 
tarded Green's functions. These are found by computing the quasinormal frequencies of the 
bulk theory [12l []]. As the characteristic oscillations of black objects, quasinormal modes 
typically decay exponentially in time (see [13] for an example and [14J for a review). By 
searching for modes whose frequencies have instead a positive imaginary part, one can study 
the linear stability of the bulk theory and thus the likelihood of phase transitions in the 
boundary theory. For example, the quasinormal frequencies for shear-type electromagnetic 
and gravitational perturbations of a Reissner-Nordstr0m AdS4 black hole were computed in 
[T5] and no such instabilites were found. Many studies, even those in which the retarded 
Green's functions themselves are computed, often restrict to the regime of long wavelengths 
and/or small frequencies. In this limit, the effective description is given by hydrodynamics 
and only the leading poles are significant [16] . Some examples are [HI [H] and [19] Q 

In this paper we aim to go beyond these studies. We compute the retarded Green's 
functions in full for conserved currents in the shear channel, considering arbitrary frequencies 

1 See also [20 for a review of the recent developments in applying holographic methods to understand 
non-equilibrium physics in strongly-coupled field theories. 
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and momenta. We consider a field theory which is assumed to be dual to an electric Reissner- 
Nordstr0m AdS4 black hole with planar horizon. On general grounds, this field theory is 
expected to live in (2+l)-dimensional Minkowski space and have a U(l) charge density. We 
will work exclusively at non-zero temperature by studying a non-extremal black hole. This 
calculation extends the work of [T5] . 

We wish to explore the detailed structure of these correlators. In this paper we show 
that their poles have varying residues and that this information can be used to assess the 
dominance of these poles beyond the hydrodynamic regime. After analysing the properties of 
these correlators for large and small values of the boundary theory parameters, we examine 
the curious and intricate motion of their poles at intermediate parameter values. 

In the following section we outline the bulk theory and the prescription for computing 
the correlators of interest. For more details we refer the reader to [T5] . whose conventions 
we share. In Section [3] we present and discuss our results and the various consistency checks 
they satisfy. We conclude with thoughts for future work. 

Note added: For comparison, we show in Appendix |C] that we recover the hydrodynamics 
of the dual theory if we take the appropriate limit of our more general numerical results. We 
also illustrate how the effective theory provided by hydrodynamics breaks down for shorter 
wavelengths and larger frequencies. Our analytical results agree with those of [21J. 

2 Background and method 

An action for a U(l) gauge field A coupled to gravity with a negative cosmological constant 
in four spacetime dimensions is 

S = 2^f / d ' X ^ { R+ J?~ L2F ^ FGb ) > ( 21 ) 

where L is the curvature radius and F = dA. A solution to this theory describing an electric 
Reissner-Nordstr0m AdS4 black hole with planar horizon can be written 

2 r 2 

ds 2 = ^ (-fdt 2 + dx 2 + dy 2 ) + ^7dr 2 , A = (jl(i-j > ) dt, (2.2) 



where 



L 2 v a ' r 2 f 

Qr 



/(r) -!-(! + <« + (2.3) 

The (outer) black hole horizon is at r = ro and has temperature 

Fluctuations about this background can be expressed in terms of gauge-invariant master 
fields $± [22]. When written in terms of these fields, the linearised equations of motion 
decouple. We are interested in the retarded Green's functions of operators in the shear 
channel of the boundary theory. After using rotation invariance in the (x, y) plane to set the 
momentum in the y direction to zero, we can focus on correlators of J y) T xy and T yt . The 
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master fields appropriate to this case, as well as further references and details, can be found 
in [15]. 

To aid numerical work, we introduce the dimensionless coordinate z = r /r. In these 
coordinates the AdS4 boundary is at z = and the black hole horizon is at z — 1. The 
equation for the master fields, given in [15], becomes 

z 2 f(f& ± )' + [-zff + Q 2 z 2 (w 2 - q 2 /) - 2Q 2 9±z 3 f] S± = 0, (2.5) 
where a prime denotes differentiation with respect to z and 



3 / 1 \ / / 16 / 1 x 2 



*W = A 1 + cF^ 1± V 1 + n 1 + <F 1 q2 '- (2 - 6) 

The dimensionless parameters vo and q are respectively the frequency and spatial momen- 
tum of the perturbation in the boundary theory, as measured with respect to the chemical 



potential. The latter is identified from the bulk theory to be /i, as given in (2.3). 

To obtain the retarded Green's functions of interest, we must adapt the original prescrip- 
tion of [1] to the case of coupled perturbations, considered in [22 EU [25] . We must first solve 



(2.5) numerically, subject to infalling boundary conditions at the horizon. Then we use the 



asymptotic fall-off of $±(z) at the AdS4 boundary, given by 



$± -°$ ± (l + n±z + ...), (2.7) 



to extract n±(tx>,q) from this solutionjj The matrix of retarded Green's functions in the 
shear channel can be constructed from these functions. Whilst the method can be found in 
[15], we present the results generalised to the non-extremal case here: 

j 2q 2 Q 2 (fl_n + -fl + n_) 
Gyt,yt - (2.8) 

_ 2w 2 Q 2 (g.fl + -g + fl.) 

3(g+-g-) { ' 

G xyyt = -^ Q2 ^- 9 ^ (2.10) 

_ q « P Q 2 (g_n + -g + n_) 
_ 2 qo >Q 2 (n + -n_) 

Gy ' xv = 2 T(1+ - gV (2 ' 13) 



2 More details and a numerical issue can be found in Appendix B 
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2g 2 (n + -n_ ) 

Ka+ - 9-) 
8 (g ± n+-^_n_ ) 



— — + (2-15) 



Note that a contact term has been suppressed in the G yt , y correlator and that all correlators 
are functions of to and q . Note also that we can pick the normalisation of two-point functions 
because we have the freedom to choose an operator basis on the boundary. As such, we have 
chosen to normalise all correlators with respect to G Vj y. These results reduce to those given 
in [15] for zero temperature, i.e. for Q 2 = 3. 

We set L = 1 in our numerics and absorb a factor of r into (£, x, y). As a result, \i — Q. 
In the rest of this paper, we refer to II ± as 'the retarded Green's functions for $±' for brevity. 

3 Results and discussion 

To ensure that the retarded Green's functions for $± have been calculated correctly, we must 
check that the locations of their poles match the quasinormal spectrum for the appropriate 
bulk fluctuations [T2l [1] . In Figure [j] we demonstrate that this is indeed the case for our re- 
sults. The quasinormal frequencies were computed using the determinant method pioneered 
in [26] and explained in detail in [27jj^] We find precise agreement with the quasinormal 
frequency plots for non-zero temperature shown in [T5] , 

First, a note on our terminology. From the gravity perspective, modes with quasinormal 
frequencies on the negative Imvo axis are associated with purely decaying perturbations of 
the black hole. We refer to these as 'on-axis' modes. The 'off-axis' modes have an additional 
oscillatory component. 

As an example, we present in Figure [2 the full result for the G xy ^ y correlator, given in 



(2.12). It can be seen from (2.8 2.16) and (2.6) that the poles in the various correlators are 



determined completely by II±, which we now study for various parameter ranges. 
3.1 Residues 

In Figure [i] we show surface plots of fl± on the complex vo plane. They verify the usual 
symmetries of retarded Green's functions, namely 



n ± (n>,q) = n ± (-n>,q), (3.1) 



which can be shown to follow from the vo <-> —vo symmetry of the equation of motion (2.5) 
and the infalling boundary conditions. 

Note that the off-axis poles have much wider peaks than those on-axis, indicating larger 
residues. We would like to understand how important these residues are in determining the 
behaviour of the boundary theory. Suppose we assume that II± have only simple poles which 



3 SAG would like to thank Sean Hartnoll for useful tips on generating these quasinormal spectra. 
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Figure 1: A comparison between density plots of |fi±| on i/ie complex vo plane (left panels) 
and the quasinormal frequencies for $± . T/ie io^> row is for $ + and i/ie bottom row is for 
All plots have q = 1 andTf/a = 0.09. we discuss later, the on-axis modes are weaker 
but have been tested thoroughly against the quasinormal spectrum in a finer plot. 



are well-separated (an assumption we question in a later subsection). Then they will have 
the form 

ILJto.q) = > — . - + analytic pieces, (3.2) 

^tx> -»i(q;T//x) 

where the locations of the poles to^ in the complex vo plane have real part rb^(q;T//i) and 
imaginary part — I^(q;T//x). By performing a Fourier transform, we show in Appendix |A| 
that the quantity 

determines the contribution of a given pole to Ii±. Beyond the hydrodynamic regime, i.e. 
for t — t f < 1, the residues can therefore be cruicial in determining which poles dominate the 
linear response. 
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Figure 2: ^4 comparison between a surface plot of \G xy ^ y \ on the complex to plane (left) and 
the appropriate quasinormal frequencies. Both plots have q = 1 and T j \i — 0.09. 

3.2 Behaviour at large frequency and momentum 

At large to and q we can perform a WKB analysis of the equation of motion. This yields 

f[ = ±iQ^/w 2 -q 2 (3.3) 

for both f[±, where we consider real to. So for to > q or q > a) the growth of \fi±\ is 
linear in the larger parameter, which we demonstrate in Figure [4} We note the zeroes at 
vo = ±q for large Revo and q in confirmation of our resul t an d also the strong presence of 



3.4 



the hydrodyanmic pole in to be discussed in a Section 

We have found the full retarded Green's functions for a non-trivial density matrix in the 
boundary theory. To isolate the modification from the vacuum, i.e. the effect of adding the 
black hole to the bulk theory, one can subtract off this asymptotic behaviour. However, 
this comes at the price of introducing the square root branch cut shown above. The cut 
originates from the retarded Green's functions of the CFT dual to pure AdS4, a fact which 
has been obscured somewhat by moving to master fields. 

3.3 Spectral function 

For a field theory in thermal equilibrium, the spectral function gives the number density of 
states in the ensemble with a particular vo and q. The spectral function is defined up to 
normalisation by 

p(t»,q) = -ImGfl(»,q) (3.4) 
for real to, where Gr is the appropriate retarded Green's function. 



The features on the Re to axis can be explained by studying the general form (3.2). 
Making use of the symmetry ( |3.1[ ) we consider a straightforward re- writing of this expression, 
with dependence on q and T/fi suppressed: 

n±(to) = > — — : — — : — +> J —— + analytic. (3.5) 

v J \xo -\xOi\+iTi xo + \xDi\+iTi\ ^w+iTj 
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Figure 3: Surface plots of U + (top row) and fl_ (bottom row) on the complex to plane at 
q = 1, T/n = 0.09. We show (from left to right) the real part, imaginary part and absolute 
value of each. All plots have the same orientation on the plane, as indicated. 



Here, i runs over one half of the off-axis poles and j runs over the on-axis poles. Taking real 
and imaginary parts of the above form we obtain 



Refl ± = 



aiW - (ai\vbi\ - biTi) ciiVt + (a;|t»;| - 6^) 



(tt>-|«M) 2 + i1 



(» + l* 4 |) 2 + r? 



+ ^2 , + analytic 



» 2 + n 

3 3 



hvo - (aiTi + bi\Wi\) hw + (aiTi + bi\Wi\) 



(w - \vbi\y + r 



(w + jwiir + r 



+E 



Cj w 



n> 2 + T2 

3 3 



(3.6) 
analytic. 
(3.7) 



In Figure [5] we have taken a slice of the (Reu>,q) plane at q = 10 -6 . The expressions above 
take their largest values whenever to = or ±|tt^|, which is indeed what we find in Figure [5j 
In effect, the presence of poles lower down in the complex to plane is 'projected' onto the 
retarded Green's functions at real to. 

Note that in Figure[5]we have subtracted the large to and q behaviour discussed previously 
in order to reveal the peaks that occur away from to = 0. Our numerics confirm that the 
modified retarded Green's functions tend to zero as Re to increases. We have chosen to 
subtract off the positive branch. 

By considering q <C 1 we are working in the long- wavelength regime. If instead we 
consider q 3> 1, for which the spatial perturbations are much smaller than the scales in the 
boundary theory, the poles on the Im to axis move much further down the complex to plane 
and their effect on the spectral function diminishes. 



Figure 4: Surface plots o/|n+| (left) and |EL| (right) on the (Reto,q) plane atT/fi — 0.09 to 
demonstrate the large vo andq scaling. Note that we avoid q =0: in that case the symmetry 
is enhanced and so the master fields are different, as discussed in [T5] . 

3.4 Pole motion 



As noted in [15], the lowest quasinormal frequency for $_ approaches zero as q — )• 0. We 
found the same behaviour for the lowest pole in II_; as such, this can be identified as the 
hydrodynamic pole. The behaviour of the poles as q and T/fi are varied is particularly rich 
and, to our knowledge, has not been studied in detail before. 

The trajectories of the on-axis poles of f[_ at fixed T/fi are shown in Figure [6j As q 
is increased, all these poles move down the Imro axis. The on-axis poles of II + have very 
simple trajectories and so we focus on the behaviour of 

Firstly we note the apparent continuity of the white line characterising the hydrodynamic 
pole. The physical implication of this behaviour is that in the hydrodynamic limit, where 
this is the only relevant pole, the dispersion relation of the quasi-particle corresponding to 
this pole is a smooth function of q. This is given by 

w = -ifiD^ 2 (3.8) 

with diffusion constant In the left of Figure [6] we can see there are values of q for which 
a given on-axis mode disappears, as indicated by a gap in the associated trajectory (for 



example at q w 1.6 and to « — 0.6i). If (3.2) holds, then we must conclude that the residue 
vanishes at these (to, q). With this in mind we see that, for finely tuned q, it may be possible 
for a previously sub-leading pole to give the dominant contribution to II_. Note that whilst 
the quasinormal plots can show this continuity, full knowledge of the correlators is required 
to determine which pole dominates. 

Also, we observe that the trajectories of all on-axis poles become less pronounced as q is 
increased, which indicates that the residues of the poles become smaller. Furthermore, there 
are no poles at small |tx>| for sufficiently large q. 

The 'interactions' between the poles of II_ are quite intricate as T j \i and q are varied. 
These features are best displayed in animations, but we give a brief summary herej^] 

See Appendix |c| for a discussion of the diffusion constant. 



5 Animations are available at http : //www . maths . dur . ac . uk/ ^rcqn58/ research . html These show the 
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Figure 5: Slices through the (Reto,q) plane of TI + ^op rou^ and II_ (bottom row) aiq = 10 6 
and T/ 11 = 0.09. T/ie /e/i and panels show the real and imaginary parts, respectively. 

We observe two types of behaviour: 'repulsion' and and 'clover-leaf crossing'. Suppose we 
fix T//LL and increase q. All poles on the Imvo axis move downwards but the hydrodynamic 
pole moves down further for the same change in q , as shown in Figure [6| It slows down as it 
approaches the next pole and the latter speeds up but the two never touch, just as for two 
like magnetic poles being brought together. As this second pole approaches the third, the 
two coalesce into a single object, split and move off- axis, coalesce a second time then split 
and move on-axis in a clover-leaf pattern. 

Now suppose we fix q at some larger value and increase T j \i. Again, all the poles move 
downwards as T j \i is increased, except there is always one pole which moves upwards. An 
example of this is shown in Figure [7] and a sample crossing is depicted in Figure [5J We also 
see several repulsions, so it appears that a different pole moves upwards each time. 

We have so far not been able to tag the poles in a clover-leaf crossing and track their 
motion independently. The fact that the poles coalesce follows from the symmetry (3.1): 



the poles must be symmetric about the Imto axis, so they cannot move round each other 
but must cross. However, we do not know if this forms a sum of simple poles or a double 



pole. If the latter occurs then (3.2) must be modified, so it would be interesting to map the 



behaviour of II_, or even just the hydrodynamic pole, over a larger region of (q, T / jj) space. 

From the gravity perspective, we have shown that certain quasinormal modes switch from 
purely decaying to partially oscillatory behaviour as the parameters are varied. For example, 



motion of the quasinormal frequencies for We observe the same behaviour for II_ and as such it is 
not an artifact of using a matrix of finite size in the determinant method. These animations also show the 
behaviour of the off-axis poles, which remains to be understood. 
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Figure 6: Trajectories on the (Imto,q) plane of the on-axis poles of II- (left) and the motion 
of the associated quasinormal frequencies. We fixT/fi — 0.09 and consider 0.1 < q < 5. 

for the central plot of Figure [8] there exist two modes degenerate in energy which move in 
opposite directions. Now, those modes closest to the Re to axis dominate the late-time 
response of the black hole. By studying this crossover we are exploring the intermediate 
response for perturbations of various q of black holes with different T j \i. From the field 
theory perspective, the dispersion relation for the excitations associated with these poles 
develops a real part at this crossover. 

4 Outlook 

In this paper we have considered a field theory which is assumed to be dual to a non- 
extremal Reissner-Nordstr0m AdS4 black hole. We have computed numerically the full 
retarded Green's functions for conserved currents in the shear channel of this theory for 
non-trivial density matrix. Our results are consistent with those of [15J, but we have gone 
beyond that treatment to show many interesting features. In particular, the linear response 
of the boundary theory exhibits rather peculiar behaviour beyond the hydrodynamic regime 
which we would like to understand in more detail. 

An immediate application of our methods is to study the sound channel correlators of this 
non-zero temperature field theory. This case was not considered in [28]. We would also like 
to compute the full retarded Green's functions for conserved currents in the shear channel 
at zero temperature. We expect qualitatively different behaviour at zero temperature. It 
has been argued in many examples (see [27], for instance) that the double zero in fir) at 
the horizon of extremal Reissner-Nordstr0m AdS4 leads to a branch cut in the appropriate 
retarded Green's function(s) of the boundary theory. This was first seen explicitly at small 
frequencies in [8]. Whilst numerical results were reported in [28j for the sound channel at 
zero temperature, only the small-to behaviour was reported in [15] for the shear channel, 
based on the analytical methods developed in [8]. We would like to go beyond this regime, 
for which numerical methods are required. However, the double zero in fir) shows up as an 
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Figure 7: Trajectories on the (Imto, T j \i) 'plane of the on-axis quasinormal frequencies of 
We fix q = 3 and consider 0.1 <T/ji< 0.25. 

irregular singular point in the master field equations, leading to difficulties in the numerical 
computation of the retarded Green's functions by our methods. It would be interesting to 
adapt our methods to tackle zero temperature for the shear channel. 

Note that we are taking a 'bottom-up' approach to the AdS/CFT correspondence. As 
such, the boundary theory is not well-defined and we may miss some interesting features, 
particularly at low temperatures. The hope is that such toy bulk theories can be embedded 
in a consistent truncation of ten- or eleven-dimensional supergravity. Important milestones 
in this regard are [29j [30J [3]] . We believe our techniques will be useful in exploring such 
setups. 
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A Retarded Green's functions in position space 



In this Appendix we derive the position space representation of the general form for II 



± 



given in (3.2). First, consider the simpler function 



ft ± (»,q) = JjQ- — ( y (A.l) 

iv - t»(q) + ?r(q) 
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Figure 8: Density plots of |n_| on £/ie complex vo plane at fixed q = 3 are shown in the upper 
row. The values ofT/fi are (from left to right) 0.185, 0.2 and 0.215. These plots have been 
cropped to include only the poles which move off- axis and the leading poles. For comparison, 
the quasinormal frequencies at the same values of parameters are given in the second row. 
The final row contains two quasinormal frequency plots, one at T/fi = 0.185 and one at 
T/fi = 0.215. These show a larger region of the complex vo plane, including the origin, for 
orientation. 
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where we have suppressed the dependence on T//j for clarity. Fourier transforming with 
respect to u and k we find 



= A* 



II : (.r ,•')-// / ^| /* (<t) . f '' ^ 

^ a6 ^ ^fa) e -i(a+w(q )-i/*r(q )) (t-t')+ifc-(x-x') 

(2tt) 2 a 

/dfc , Trfa f 



27T a 



_ e -ia{t-t') 



We can integrate the factor in square brackets using contour integration. First we nudge the 
pole a short distance into the lower half-plane. The contribution along the curved section 
of a semi-circular contour will vanish for different ranges of t — t' . We recover the following 
result for the above integral over a: 

lim ! ^^— e -Ht-t') = lim { -27nlinw ie ie^), (t - t> ) > 

e^O J 27T a + it e^O \ 0, (t — t') < 



lim 



-ie-(*-*') £ , (t - t') > 

"-io \ o, (t — f ) < o 

i6(t-t'). 



U ± {x - x') = -iyB{t - t') I _^i ? ( q ) e -Mr(c,)(t-*') e --(c.)(*-*')+*(^'). ( A .2) 



Thus, 



The retarded Green's function of ( |3.2[ ) is a sum over poles. We apply the above procedure 
to each term in that sum to obtain 

fl ± ( x - x >) = - ifI 6(t -t>) /(gjE [^(q)^^^^^'^-^^^*-*')] e^*-*). (A.3) 

It is clear from this result that the contribution of e ach pole to fl± is characterised by the 
quantity R i (c\)e~ fiTi ^^ t ~ t \ as claimed in Section 



3.1 



B Numerical method 

In this Appendix we give details of the method used to extract the retarded Green's functions 
for $± numerically. As outlined in Section [2| we must first solve (2.5). This is a second-order, 



linear ODE with a regular singular point at z = 1 (the horizon). We choose the following 
ansatz, where $ denotes one of the $±: 

$0) = (z-l)- M ^ T (f)(z). (B.l) 

The first factor imposes the infalling boundary condition at z = 1. A unique solution is then 
specified by 0(1), which we are free to choose because the equation is linear. We expand 
(j) about z = 1 up to some order, AT, to generate the initial condition for a Runge-Kutta 
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algorithm at some z — 1 — e. To extract n±(ro,q), we integrate out to the boundary and 



match this numerical solution to the asymptotic expansion given in (2.7). The matching can 
be performed using a root-finding method, for example. 

The main numerical issue comes from the initial condition. Naively, we would like to 
choose a small e so that the series expansion of (j) is accurate with only a few terms. However, 
we found that the pattern of poles was completely washed out below a certain line in the 
complex vo plane. This numerical instability appears if e is small for Imto large and negative 



(and/or for T / ji small). To see why, note that the ansatz (B.l) for $ at z — 1 — e contains 
the factor 

This factor becomes very large in these regimes if e is small, leading to round-off errors. 

Thus, a larger e must be chosen in these regimes. As a consequence, a sufficiently large 
TV must be chosen to offset the error from starting the integration further from z — lj^] 
The values of e, TV can be constrained by ensuring that the locations of the poles match the 
quasinormal spectrum for the appropriate bulk fluctuations, as stated in Section (3) 

C Hydrodynamics 

In this Appendix we show that we recover the correct hydrodynamics if we take the appro- 
priate limit of our numerical results. We begin by reproducing the analytical expressions for 
II± found in [21] for the benefit of the reader. 

As discussed in Section^ the matrix of corrrelators can be constructed from fl±. We 
now compute the latter to leading order in vo and q 2 in the spirit of [T7j Q Firstly, as in [21], 
we consider the following field redefinition, where $ denotes one of the $±: 

$(z) = (z- l)- lV0 ^ T F(z)G(z). (C.l) 

The aim is to choose G such that F tends to an vo- and q-independent constant, F , in the 
hydrodynamic limit ro,q 0. This choice greatly simplifies the subsequent analysis, as we 
will soon show, and leads to the equation 

zfG" + zf'G' - (/' + 2Q 2 g ± (0)z 2 )G = (C.2) 

which has solutions of the form 

for*. 



G(z) = a(z + b), where b=< 4( ? 2 + ' (C.3) 

I for <I>_. 

Here, a is a constant we can specify by a choice of FG at the horizon. The resulting equation 
for F is 



F +i—i + j + ir ]F + 



a fa-l f 2G' 



+^K-q7)-2Q^ ± (q, / 



z-l\z-l f G 

= 0, (C.4) 



6 SAG would like to thank Patrick Kerner for discussions on this issue. See [32] for further evidence. 
7 There can be no odd powers of q due to parity. 
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where we have defined 

aiXO) = ~Br and 9±to)=9±{0-9±$)- (C5) 

Now that we have defined G in this way, we can write F as a series in vo and q 2 , i.e. 

F(z) = F + xoF^z) + q 2 F 2 (z) + 0(n> 2 , n>q 2 ), (C.6) 

then collect terms in ( C.4[ ) at each order. Before moving on to find F\ and F 2 , we note two 
points which guide the analysis. First, the inf ailing boundary condition at z — 1 implies 
that F be regular there, and second, we only require the small z expansion of F in order to 
read off tl± (as described in Appendix |b| . 

Collecting terms at 0(tx>), the equation for F\ is 

f; + (7 + it) f; + ^T)(— + 7 + ir) =0 ' <a7) 

which can be integrated once to give 

^' = -— P-T + ^Y (C8) 



vo \z-l fG 2 

The constant ci is determined by demanding that the 0(z — terms cancel, and is found 
to be 



(3-gj) 



2\3 



ci = < i6g4 + ' (c.9) 

^(3-Q 2 ) for$_. V ; 



for $ 



The factor of a 2 has been extracted because F cannot depend on a. 
Collecting terms at 0(q 2 ), the equation for F 2 is 



n +{j + 2 -§-)n- ^(1 + 2~g ± (l)) = 0, (CIO) 



which can be integrated once to give 



2 f(z + b)* 



c 2 + b[bz + z 2 + 2g±(l) I — — I — — + — + 



bz 2 2z 3 \\ z 3 g±(l)z^ 



2 3/73 2 



(C.ll) 



The constant c 2 is determined by demanding that the 0(z — 1) 1 terms cancel, and is found 
to be 

( 27+63Q 2 +29Q 4 +9Q 6 r ^ 

"* ^ 1 for$_. V ; 



3(1+Q 2 ) 



To extract we need only keep terms up to 0(z) in $±. Let's focus on II + for now. 
Expanding each factor to this order we obtain 



$ + = (-1) Q (1 - az)a [F + 0)^(0) + F[(0),z) + q 2 (F 2 (0) + ^(0)*)] (z + b) 
+ 0(z 2 ,t» 2 ,t»q 2 ). 
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By comparison with (2.7) and substituting for b we obtain 

iF[(0) 



n. 



4Q 2 

'3(1 + g 2 ; 



+ IVO 



Q 
3-Q 2 



F 



+ q 2^M + 0(tt> 2 ,tt>q 2 ). 



(C.14) 



Substituting for F-[(0) and ^(0) we obtain the final result 



n, 



4Q 2 Q(3-Q 2 ) 2 2 Q 2 (27 + 63Q 2 + 29Q 4 + 9Q 6 ) „. 2 2 . , 1 . 
+M»-^ LA : zr^r-M-^ r^r: zttI -^+0(t» , t»q 2 )- (C.15) 



3(1 + Q 2 ) 9(1 + Q 2 ) S 



27 {I + Q 2 Y 



The case of Il_ is slightly different because the ^-dependence of F\ and F 2 is different. 
We find 



Fi 



aF fQ 2 -3 



+ d 1 )+0(z), 



F 2 = Q 2 F 



1 



,3(1 + Q 2 )z 

where d\ and d 2 are constants. Using b = we obtain 



+ d 2 +0{z), 



(C.16) 
(C.17) 



$_ = (-l) Q aF 
and thus 



a(3 - Q 2 ) 



n_ 



q 2 g 2 



3(1 + g 2 ) 

1 — ac?i + <\ 2 Q 2 d 2 
a(3-Q 2 ) + 4 



+ (1 - ctdi + <\ 2 Q 2 d 2 )z 



+ 0(z z ,w z ,w^) (C.18) 



+ 0(t» 2 ,t»q 2 ). 



(C.19) 



3(l+<3 2 ) 



We now provide a check on our hydrodynamic result (C.19). In Section 3.4 we identified 
the lowest pole of II_ as belonging to the hydrodynamic mode. The location of the pole gives 
the dispersion relation for this mode. Fluctuations transverse to the direction of momentum 
flow, as in the shear channel we are studying, excite diffusive modes. As such, these have a 
dispersion relation of the form 

to = -i^Dq 2 (C.20) 



with diffusion constant D. This is indeed what we find in (C.19) if we identify 

D * 



3(1 + Q 2 )r ' 



(C.21) 



Now, in a hydrodynamic system which is also conformal, we expect the following relations 
to hold for energy-momentum transport: 



?7 e 
D = — with P=-, 



(C.22) 



where 77, e, P are respectively the shear viscosity, energy density and pressure of the system. 
Using the thermodynamics of our black hole (as quoted in [15]). 



K 2 L4 



(C.23) 
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we obtain 



V 



1 

47r' 



(C.24) 



This result is expected on general grounds for a theory dual to a two-derivative gravity 
theory [33J. 

Note that this D can only be interpreted as the diffusion constant in the regime where 
hydrodynamics gives an effective description of the dynamics, i.e. for T > /i. Outside this 
regime, we expect a different effective theory to be valid. In particular, we should not trust 
(C.21) for T <C \i\ either the dispersion relation needs modifying, or the other poles in II_ 
cannot be ignored, or both. 

We can calculate the diffusion constant numer ically by fitting the location of the lowest 
quasinormal frequency of II _ to the form (C.20). Our numerical results agree with the 
analytical result in (C.21) for T 3> /x, as shown in Figure [9| To reiterate: to calculate D 



we only need knowledge of the lowest quasinormal frequency and not our more complete 
numerical results for the correlators. 




Figure 9: Plot of jiD as a function of T/ji. The solid curve is the analytical result from 



(C.21) and the points were extracted from the motion of the lowest quasinormal frequency. 



Now we show how this agreement breaks down beyond the hydrodynamic regime, and 



as such our results are essential. In Figure 



10 



we get the correct contact term for II + as in 
( |C.15D and the results are in good agreement over a small range in to, but they soon deviate. 
The agreement persists over a larger range of vo if we work at a larger value of T//ij^] 
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